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Previously developed analytic models for the evolution of cosmic string and monopole networks are 
applied to networks of monopoles attached to two or more strings; the former case is usually known as 
cosmic necklaces. These networks are a common consequence of models with extra dimensions such 
as brane inflation. Our quantitative analysis agrees with (and extends) previous simpler estimates, 
but we will also highlight some differences. A linear scaling solution is usually the attractor solution 
for both the radiation and matter-dominated epochs, but other scaling laws can also exist, depending 
on the universe's expansion rate and the network's energy loss mechanisms. 



I. INTRODUCTION 

Currently favored fundamental theories suggest that 
we live on a four-dimensional slice of a higher- 
dimensional universe, and while most forces are confined 
to our slice, gravity may leak off. These so-called brane 
world models may provide natural explanations for in- 
flation (here known as brane inflation), and topological 
defects may at the end of it. The type of defect network 
which is formed and its basic properties (from its energy 
scale to whether or not it is long-lived) will depend on 
the specific details of the model in question [l| . A thor- 
ough overview of the subject may be found in the book 
by Vilenkin and Shellard Q and in more recent review 
articles 

Most of the past work on defects concerns the simplest 
models of cosmic strings, on the grounds that they are 
cosmologically benign, and are a generic prediction of in- 
flationary models based on Grand Unified Theories [(| 0] 
or branes [8, 9], while domain walls and monopoles tend 
to be cosmologically dangerous and tightly constrained. 
However, it is clear, particularly in the context of models 
with extra dimensions, that networks containing more 
than one type of defect will often be produced. Two 
examples that have attracted considerable interest are 
semilocal strings 

CHI and cosmic necklaces [H}. For 
the latter these claims have been made both in the con- 
text of brane inflation and in string theory itself 

This is the third report on an ongoing project which 
is addressing some of these issues. In the past we have 
developed [19| an analytic model for the evolution of net- 
works of local and global monopoles (23413 ■ The model 
is analogous to the velocity- dependent one-scale model 
for cosmic strings [23l - [25| . which has been extensively 
tested against field theory [2l| H3] and Goto-Nambu sim- 
ulations [13, H3 ■ This was then extended [H| to the case 
of monopoles attached to one stri ng ( the so-called hybrid 
networks [13]), as well to vortons [3l|,[32|]. Here we study 
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defect networks where monopoles are attached to two or 
more strings. Networks of the first type are commonly 
called cosmic necklaces; we shall refer to networks of of 
monopoles attached to three or more strings as cosmic 
lattices. 

The behavior of necklaces and lattices is qualitatively 
similar, and will be for the most part treated together in 
this paper, through we will point out the small existing 
differences. However, their evolution differs in several key 
aspects from both that of individual monopoles and that 
of monopoles attached to a single string (usually called 
hybrid networks) [III Hsj. The main difference is that 
necklaces and lattices form stable, long-lived networks 
which usually reach a scaling solution. 



II. COSMIC NECKLACES AND LATTICES 

We start with a brief overview of previous results on 
the evolution of necklaces and lattices. This is by no 
means exhaustive; the aim is to highlight the dynamical 
aspects we will need to model. A more detailed discussion 
can be found in [lj] as well as in other earlier references 
that we will point out where appropriate. 

The defect networks of interest form via the symme- 
try breaking pattern G — > K x U(l) — > K x Z^. If 
G is a semi-simple group, the first phase transition pro- 
duces monopoles while in the second each monopole be- 
comes attached to N strings. If K is trivial all the 
(Abelian) magnetic flux of the monopoles is confined into 
the strings, and there are no unconfined fluxes. How- 
ever, unconfined non-Abelian magnetic fluxes can exist 
in the generic case. The previously discussed hybrid case 
[H| corresponds to N = 1. Here we will discuss the 
case N > 2; N — 2 corresponds to cosmic necklaces and 
N > 3 to cosmic lattices. 

The corresponding defect masses will be m ~ (Ait je)r\ m 
and \i ~ 27r?7g, while the characteristic monopole radius 
and string thickness are S m ~ (ery. m ) _1 and S s ~ (er/ s ) _1 . 
There are also scenarios where the intermediate phase 
transition is absent, G — > K X Zn, in which case an anal- 
ogous network still forms but the role of the monopoles 
is now played by solitons that are usually called 'beads'. 
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In this case the two energy scales are obviously similar, 
that is r) s ~ r) m . 

Up to the second transition (if it exists) the models for 
plain monopoles [l^| apply, but once the strings form a 
separate treatment is needed. Two key differences are im- 
mediately apparent. The evolution of isolated monopoles 
can be divided into a 'free' (pre-capture) and a post- 
capture period, with captured monopoles effectively de- 
coupling from the network and losing energy radiatively 
until they decay (this is analogous to the evolution of cos- 
mic string loops). In the present context the monopoles 
are effectively captured by the strings, and one needs to 
account explicitly not only for radiative losses (for ex- 
ample gauge radiation if there are unconfined magnetic 
charges) but also for the force the strings exert on the 
monopoles — depending on the context, the forces due to 
the string(s) or the other monopoles may be the domi- 
nant ones. 

If all the strings attached to each monopole have the 
same tension (which we will assume to be the case in 
the present paper) then all the strings pull it with equal 
forces, and therefore there is no tendency for a monopole 
to be captured by the nearest antimonopole, unless their 
separation is of order 5 S . If there are N strings attached 
to each monopole, its proper acceleration is given by the 
vector sum of the tension forces exerted by the strings. At 
a back-of-the-envelope level, each force is of order / ~ /i, 
and hence one expects that a ~ \ijm. Monopoles should 
therefore be accelerated to relativistic speeds provided 
that the characteristic length of string segments, L s , is 
such that pL s 3> m, that is L s /S s 3> T) m /r) s . 

Aryal et al. |33j ] first studied the formation and statis- 
tical properties of these networks, for N = 2 and N = 3, 
showing that for N > 3 a single network is formed. In 
all cases they find that the system is dominated by one 
infinite network comprising more than 90% of the string 
length. Some finite networks and closed loops do exist, 
in numbers rapidly decreasing with their size. Finally, 
most of the string segments have a length comparable to 
the typical distance between monopoles (much larger seg- 
ments being exponentially suppressed). This justifies our 
assumption of an inter- monopole separation, L m , compa- 
rable to L s . 

The cosmological evolution of these networks was first 
discussed by Vachaspati & Vilenkin [34[ , who argued that 
assuming that the radiation of gauge quanta is the dom- 
inant energy loss mechanism of the networks, they reach 
scaling with a characteristic lengthscale L ~ (Vs/Vm)^ 
and the monopoles become highly relativistic. These net- 
works can also lose energy by producing closed loops of 
string and small nets. The effect of these is harder to 
estimate, but as we shall discuss it is fairly easy to model 
phenomenologically. . 

Specifically they divide the network's energy into 
string and monopole parts, with p m = ftp, p s = (1 — j3)p 
and an effective equation of state 3p/p = j3+(l — ft)(2v 2 s — 
1) with V s being the string velocity and assuming v m ~ 1. 



The evolution equation is then 

p = —3H(p + p) — nw (1) 

where w ~ (ga) 2 /6ir for gauge radiation losses and 
L ~ n -1 / 3 is both the length of string segments and 
the average monopole distance. Assuming a self-similar 
evolution one can also set p s ~ /in 2 / 3 . 

For the case without unconfined magnetic fluxes (that 
is, no radiation), they claim L oc ( Q , expecting a < 1: 
there is no scaling and the defects eventually dominate 
the energy density of the universe. But note that in 
saying this they are specifically thinking of the radia- 
tion epoch (it is clear that the behavior of their solution 
depends on several parameters including the expansion 
rate). 

The specific case of cosmic necklaces has subsequently 
been studied by Brezinsky & Vilenkin [35( . They assume 
no unconfined magnetic fluxes (hence no Coulomb forces 
between the monopoles) and characterize the networks by 
a dimensionless ratio r — mj (pL) with the average mass 
per unit length of the necklaces being (r + l)p. They also 
neglect the effect of annihilations (though the validity of 
this assumption has been challenged [H, H3]), and find 
that the system tends to evolve towards large r. Again 
the necklaces are expected to evolve in a scaling regime, 
with a characteristic network lengthscale £. The force 
per unit length of string is / ~ and the acceleration 
is a -1 ~ (r + 1)£ so we expect that 

£~^rW w ~^tW (2) 

s/l + r \/l + r 

In the limit r < < 1 the monopoles are sub-dominant and 
the strings will behave approximately as ordinary ones. 
In the limit r » 1 the strings are very slow and their 
separation is small. This is a very simple toy model, as 
in fact T is generically not a constant parameter, since 
d ~ £, so these solutions are only approximate. Nev- 
ertheless, this approach has the advantage of algebraic 
simplicity, and we shall show below that in appropriate 
circumnstances it can be related to more robust mod- 
els. Analogous results have been found with a somewhat 
different toy model [171 ] . 

III. QUANTITATIVE EVOLUTION 

It's easy to start modeling these networks by u sing the 
evolution equation derived in our previous work [I9tl29j . 
Most authors at this point focus on the evolution of the 
strings, treating the monopoles (as it were) as a small 
correction. Our approach, justified in [29j], is precisely 
the opposite — we focus on the evolution of the separation 
between monopoles. 

For this context the evolution equations for the char- 
acteristic separation L and root-mean squared velocity v 
of the monopoles are 

3^ = (3 + v 2 )HL + (3) 
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dv 



Vrn 



(4) 



We have neglected the term due to friction in both 
equations (it's easy to show this is subdominant at late 
times). The energy loss term Q± in the lengthscale equa- 
tion is a renormalized quantity, accounting for the vari- 
ous losses present (including string intercommutings, ra- 
diation and annihilations), as discussed in our previous 
work. There may be a velocity-dependence of some of 
these contributions, but as we shall sec monopoles will 
typically have ultra-relativistic velocities v ~ 1 and there- 
fore this dependence can be neglected. The velocity equa- 
tion includes the force due to the strings (with a phe- 
nomenological curvature parameter k s that is discussed 
in [HI) but we have neglected that due to monopoles, 
since if it exists (which is only the case for unconfincd 
fluxes) it's always smaller than that due to the strings. 
Indeed, using the above definitions of mass scales and 
thicknesses one finds that the ratio of the two forces is 



fm km / Os i | 



(5) 



we expect the ki to be (dimensionless) coefficients of or- 
der unity, though note that they should be different for 
necklaces and lattices. 

From the velocity equation we immediately confirm 
that the monopole velocities will be driven towards unity, 
v 1, as previously stated. As for the monopole lenght 
scale, assuming a generic expansion rate a oc i A , we find 
two different regimes for slow and fast expansion rates 



L 



Q* 
3-4A 



t, A<3/4 



L cx a 4 / 3 oc i 4A / 3 , A>3/4. 



(6) 
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The former explicitly requires a non-zero energy loss 
rate — we will return to this point later. We therefore 
have linear scaling both in the radiation and matter eras 
(as generically claimed by previous authors, based on 
simpler qualitative arguments). For monopoles L a t 
corresponds to the monopole density decreasing relative 
to that of the background. However, for fast expansion 
rates the growth is superluminal, and the network will 
eventually disappear. An analogous scaling solution was 
already discussed in our previous work on hybrid net- 
works. 

Its easy to establish a link between this analysis and 
that of Vachaspati & Vilenkin [34j (which is embodied in 
Eq. [1]) and thus to carry out a more detailed analysis of 
of the possible scaling solutions. Let's apply Eq. Q]to the 
monopoles. We use both the notation and definitions of 
our previous work and those of [33 |. Since p m = ran = 
m/L 3 and the monopole equation of state is 3p = 
we get by substitution 



3— - (3 + v 2 )HL + — w. 
at r) m 



(8) 



Now, Vachaspati & Vilenkin are assuming energy losses 
through gauge radiation; noting that 



w 



(go) 2 
Gtt 



Q 



gauge 



(9) 
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we see that this evolution equation for L is exactly the 
same as Eq. [3J matching the Q terms (which in our case 
can phenomenologically account for further energy loss 
channels). 

We can also apply Eq. Q] to the strings. In this case 
2 1 3 = fi/L 2 and the string equation of state is 
l)p; again we find 



Ps 

3p 



[in 
(2oJ 



2§ -2BW + 4) 



W/J, . 



In this case we have 



(11) 



(12) 



Note the interesting fact that the dimensionless parame- 
ter Q determines the energy loss term for both the strings 
and the monopoles. The above is the usual evolution 
equation for the cosmic string correlation length (23T - [25j , 
if one assumes a constant string velocity — otherwise the 
energy loss term should depend linearly on velocity. 

The scaling solution for the monopoles has already 
been discussed. For the case of the strings the solution 
is also the expected linear scaling 



L 



2A(l + «») 



t, 



(13) 



for constant velocities and provided A(l + v 2 .) < 1. But 
this solution is an attractor, as in the case of normal 
strings: for very large lengthscales the string velocity 
would no longer be a constant (the string velocity evo- 
lution equation would drive it to smaller values), and a 
new equilibrium value with a smaller lengthscale would 
be reached. 

We can also consider more generic scaling solutions 
of the above equations, ie allowing for the possibility of 
zero energy losses (Q = 0). We will confirm that scaling 
(L oc t) generically requires Q ^ 0. Starting again with 
the monopoles, for Q ^ we have the two branches of 
the solution discussed above. For Q = the solution is 
always 



Loca 4 / 3 oci 4A / 3 ; 



(14) 



note that for A < 3/4 the lengthscale grows subluminally 
while for A > 3/4 it grows superluminally. In the absence 
of radiative energy losses, only a fast enough expansion 
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can dilute the network. We can also compare the evolu- 
tion of the monopole and background densities 

Pm Vrn t /..i 

2~Ta- ( 15 > 

For the linear scaling solution L oc t this has the form 

Pm r 1 / T?m \ 2 f T \ 2 ^ 1 
Pb Q* \mpi) \mpi) t ' 

From this we see that if gauge radiation is present then 
the energy density of the network is smaller than the 
background density. However, if the only radiative chan- 
nel available is gravitational radiation (in which case, as 
discussed in [2{J, Q+ = Q gra v ~ {Vs/ m Pi) 2 ) then the net- 
work energy density is in fact the dominant one. For the 
non-scaling branch the density is 

^ oc t 2 ~ iX , (17) 
Pb 

and the behavior depends on the cosmological epoch. No- 
tice that during the radiation era the density is a constant 
fraction of that of the background. 

These results can now be related with the toy-model 
analysis of Berezinsky & Vilcnkin. In the linear scaling 
regime we have r — p m /ps °c t^ 1 and therefore r — > 
and v — > 1 in agreement with our quantitative analy- 
sis. In the non-scaling branch r oc £ 2 ~ 4A , whose behavior 
again depends on the cosmological epoch. The radiation 
epoch corresponds to the interesting case r = const., 
while faster expansion rates (say, the matter-dominated 
epoch) dilute the monopole density relative to that of 
strings. Slower expansion rates (A < 1/2), which some- 
times occur in the very early universe in string cosmology 
models, would have a growing density ratio r. 



IV. SUMMARY AND OUTLOOK 



We have extended a recently developed analytic model 
for the evolution of monopole [l9| and hybrid networks 
[29l | to the case of monopoles attached to several strings. 
We discussed their possible scaling solutions, generically 
confirming the expectation that the network will reach 
linear scaling (with its characteristic lengthscale L oc t), 
but also showing that other scaling behaviors can occur, 
depending on the expansion rate of the universe and on 
the energy loss mechanisms available to the network. 

This completes the basic structure of analytic tools 
needed to study defect networks containing both strings 
and monopoles. These models have the advantage of 
conceptual simplicity, in addition to that of allowing a 
quantitative description of the network's evolution. Nev- 
ertheless, a further generalization will be required in or- 
der to more accurately describe the possible effects of a 
hierarchy of string tensions [17, 38, 39]. A further inter- 
esting case, which requires additional dynamics is that 
of semilocal strings [lfl, [H| > which we will address in a 
subsequent publication. 
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